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Periodic Anderson model with degenerate orbitals: 
linearized dynamical mean field theory approach 
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We investigate a multi-orbital extension of the periodic Anderson model with particular em- 
phasis on electron correlations including orbital fluctuations. By means of a linearized version of 
the dynamical mean-fleld theory, we compute the renormalization factor, the density of states, 
the spectral gap and the local correlation functions for a given set of the intra- and inter-orbital 
Coulomb interactions as well as the Hund coupling. It is found that when a certain condition 
is met for the intra- and inter-orbital interactions for / electrons, orbital fluctuations are en- 
hanced, thereby enlarging the Kondo insulating gap. This effect is suppressed in the presence 
of the Hund coupling. We also clarify how the Kondo insulator is continuously changed to the 
Mott insulator when electron correlations among conduction electrons are increased. 
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§1. Introduction 

Since the discovery of heavy-fermion materials with 
rare-earth or actinide elements, comprehensive under- 
standing of this class of correlated electron systems has 
received considerable attention. In these compounds, 
strong correlations among / electrons play an impor- 
tant role to form the heavy quasi-particle state at low 
temperatures. In particular, if correlated /-electron sys- 
tems are in the insulating phase, they are referred to 
as the Kondo insulator, which possesses a wide variety 
of compounds. To discuss the electronic properties 
theoretically, the periodic Anderson model (PAM) has 
been investigated extensively. This model is simplified 
to extract the essence of heavy fermions, which is usu- 
ally described by free conduction electrons coupled to 
highly correlated single-orbital / electrons. 

In some heavy-fermion compounds, the other pa- 
rameters dropped in the original PAM, such as the 
orbital degeneracy and the interactions among con- 
duction electrons, are also important to explain the 
experimental findings. For example, the compound 
Nd2-a;Ce2,Cu04^^"^^^ shows the unusual behavior in the 
specific heat, which may be explained by properly taking 
into account the correlations due to not only / electrons 
but also conduction electrons. Furthermore, it is claimed 
that the orbital degrees of freedom affect the Kondo- 
insulating gap around Fermi surface for another pro- 
totype of heavy fermion compound, YbBei2.^"'^'' These 
findings naturally encourage us to explore the effects be- 
yond the simple PAM systematically, such as correlations 
due to conduction electrons, the orbital effects, etc. 

In this paper, we investigate the doubly degenerate 
periodic Anderson model to discuss the effects on the 
Kondo insulator due to orbital degeneracy together with 
electron correlations for / and conduction bands. By 
exploiting a linearized version of dynamical mean field 
theory (DMFT), we show that the interplay of several 



interactions yields interesting effects on the formation of 
the Kondo insulator, which do not appear in the simple 
PAM. 

This paper is organized as follows. In §2, we introduce 
the model Hamiltonian with orbital degeneracy, and then 
briefly explain a linearized version of DMFT. In §3, we 
discuss electron correlations due to conduction bands by 
employing the single-orbital model. We then explore in 
§4 how the interplay of several distinct interactions in the 
degenerate model produces nontrivial effects on the low- 
energy electronic properties. A brief summary is given 
in the last section. 

§2. Model Hamiltonian and Method 

2.1 Two- orbital periodic Anderson model 

We consider the periodic Anderson model with two- 
fold degenerate orbitals, which may be described by the 
following Hamiltonian, 

'HpAM =^0 + Hi (1) 

<i.j> 2,m,cr 

+ ^ E (/i.c™- + ^-c-) (2) 

i.7n i.ni 

- j/J^sfi-s^- JcESa-s?2 (3) 

i i 

where fima{cima-) annihilates an / electron (conduction 
electron) with spin cr(=t,i) and orbital m(= 1,2) at 
the zth site, n^^^y = c^^^Cima and n^^^^ — fj^^^jfima- 
Here, Uj represents the hopping integral, e/ the en- 
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ergy level of the / state, and V the hybridization be- 
tween the conduction and / states. Uf{Uc) is the intra- 
orbital Coulomb interaction, C/j(C/^) is the inter-orbital 
Coulomb interaction, and Jf{Jc) is Hund coupling for 
/ electrons (conduction electrons). The spin operators 
are defined by = i Y.af3 fLa^e^pfimf} and = 
5 Z^a/3 4ma'^a/3CiTO/3, where T is the PauU matrix. Note 
that the model has two orbitals both for conduction and 
/ electrons, which are specified by the same index, m. 
This scheme has been sometimes used for investigating 
the orbital effects on heavy fcrmion systems, which may 
capture some essential properties due to degenerate or- 
bitals. In this paper, we discuss how electron corre- 
lations are developed in the Kondo insulator with multi- 
orbitals, by changing the parameters {Uf,Uj,Jf} and 
{Uc,U^,Jc} systematically. 

2.2 Linearized dynamical mean-field theory 

We make use of DMFT, which has been developed by 
several groups^^^^^^ and has successfully been applied to 
the single-band Hubbard model, the two-band 

Hubbard model,^^. 33-45) ^^e PAM^^- ^^-50) ^ -p^is 

treatment can take into account local electron correla- 
tions precisely, and thereby is exact in infinite dimen- 
sions. Then the lattice model can be mapped onto an 
effective impurity model, which has been solved by a va- 
riety of methods such as the iterated perturbation the- 
Qj.y 25,28,44) ^YiQ non-crossing approximation,^^' 3^' the 
projective self-consistent method. ^^''^'^^ Since these meth- 
ods are not efficient enough to treat the systems with de- 
generate orbitals, powerful numerical methods have been 
employed, e.g. the exact diagonalization,3'^'39.42,45) ^j^g 
quantum Monte Carlo simulation.33^35, 40,43,44) j^g^ently 
a new method to solve the effective impurity model has 
been proposed by Potthoff;'^^^ which we will refer to the 
linearized DMFT in this paper. This approach simplifies 
the procedure of DMFT by linearizing the self-consistent 
equations in the low-energy region, but still keeps the 
essential features of electron correlations. In this ap- 
proximation, the effective bath is represented by a few 
sites. In spite of this simplification, electronic properties 
in the low-energy region around the Fermi surface can 
be described rather well.^^'^^^ In fact, the critical values 
of the Hubbard model with and without degenerate or- 
bitals are in good agreement with the other numerical 
techniques. 

In this paper, we exploit the linearized DMFT, which 
may convert the original model eq. (1) to the effective 
impurity Anderson model. 



Hi, 



Ima 



+ ^/ E^ fmcrfrr 



h.c.) 



Ima 

2 2 
+ f^/ E '^mT^mi + E "mT"mi 



a, (3 a,l3 



(4) 



where Cc and Sf are the renormalized energy levels, which 
should be determined by the number of conduction and 
/ electrons in the system. Here, dj^^{dim<7) creates (an- 
nihilates) an electron in the effective bath with two sites 
(^ = 1,2). 

To determine the effective energy e; and the hybridiza- 
tion Vi , we linearize the self-energies in the small u) region 

as. 



^ff{uj) = a + buj + 0{uj^) 
E^^u) = c + dLO + olu;^) 
T.f'^iuj) = g + huj + 0{u'^), 



(5) 



where a. h.c. d.g and h are real numbers to be deter- 
mined. The renormalization factors for / and conduction 
electrons are respectively given as. 



(6) 
(7) 



On the other hand, the Green function for the lattice 
system is given as, 



(8) 



G{z) = j dkG{k, z) 

G {k,z)-\^ _y_s/c(^) z-ef-^ff{z) j- 

(9) 

We consider here the D— dimensional Bethe lattice with 
the hopping integral t = t*/VD, which results in the 
density of states, 

/'o(^) = ^V4t*2-z2. (10) 

By comparing the local Green function with the impurity 
Green function, the self-consistent equation for DMFT 
now reads 

[GoLMcc^^- (t) ^^^^ 

where W{= At*) is the band width. 

By substituting the equation (5) to the self-consistent 
equations, we end up with the hybridization and the en- 
ergy of / level for the effective bath, 



;^2 



2 



(12) 
(13) 



where M2°'' = T,j^itfj = J dx x'^po{x) is the variance of 
the noninteracting density of states and Vg = V+gis the 
effective hybridization. We note here that the quantities 
a and c in (5) can be incorporated in if and ic, and 
assume h = since ft ~ always holds in the parameter 
regime we are now interested in. 

In the following, we take = 1 as unit of the energy 
and fix the hybridization V = 1.0. We deal with the 
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symmetric case in the PAM (half-fil led bands) by s etting 
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Vi = V2 
-Uc/2 - U^, and if 



— > £1 



= 

-Uf/2- 



2 ;r 



U'f, for simplicity. 



§3. Single-Orbital Periodic Anderson Model 

We begin with the single-orbital PAM. As is well 
known, in the conventional PAM with Uc = 0, the in- 
crease of U J enhances correlations among / electrons, 
resulting in the heavy quasi-particle state. For half-filled 
bands considered here, the system naturally leads to the 
Kondo insulator with a renormalized spectral gap. Here, 
we focus on the role played by electron correlations due 
to the conduction band. To this end, we calculate the 
renormalization factors Zc and Zf as a function of the 
c-c Coulomb interaction Uc- As shown in Fig. 1 (a), the 
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Fig. 2. DOS of / electrons for the single-band PAM. 




Fig. 1. (a) Renormalization factors as a function of the c-c 
Coulomb interaction Uc- solid (open) symbols represent Zc {Zf). 
We show the results for U j = 0.0, 5.0 and 10.0 as circles, trian- 
gles and squares, respectively, (b) Local spin correlation func- 
tion (s{ ■ S?), and double-occupation probabilities {n^^n^^) and 
(<<>■ 



introduction of Uc reduces the renormalization factor Zc, 
while it does not alter the factor Zf so much. Namely, 
the c-c interaction Uc mainly renormalizes conduction 
electrons, as naively expected. We also show in Fig. 1 (b) 
the c-/spin correlation function, (Sf-S?), and the double- 
occupation probabilities, and (n|n|), for / and 
conduction electrons. Introducing U^ the double occu- 
pation probabilities of the two bands are decreased, and 
accordingly the c-/ spin correlation function approaches 
-3/4, implying that the spin singlet is formed between 
the / electron and the conduction electron. Therefore, 
the spin sector for large Uf and Uc is stabilized by the 
formation of the local singlet at each site. 

Although the interaction Uc enhances electron corre- 
lations, it may not necessarily reduce the spectral gap 
characteristic of the Kondo insulating phase. To see this 
clearly, we compute the density of states (DOS) for the 
one-particle spectrum shown in Fig. 2. For the non- 
interacting case [Uf = Uc = 0), the gap in the vicinity of 
the Fermi level is caused by the bare c-f hybridization. 
When the /-/ Coulomb interaction Uf is increased by 
keeping Uc = 0, f electrons are renormalized, resulting 
in the formation of the small gap typical for the Kondo 
insulator. In addition to the low-energy Kondo peaks, 
the large spectral weight appears around i / and ej + Uf, 
which is represented by a bunch of spiky peaks since the 
effective heat bath is represented by a few sites in the 



linearized DMFT. Similar behavior to reduce the gap is 
observed when only the c-c Coulomb interaction Uc is 
introduced, as should be expected. 

On the; other hand, when Uc is turned on together with 
Uf, the gap is increased as shown in Fig. 2. This is con- 
trasted to the above-mentioned cases possessing either 
of Uf or Uc- To make this behavior more explicit, we 
plot the gap A in Fig. 3, from which we can indeed see 
the above mentioned properties. The gap is expressed in 




Fig. 3. Spectral gap A for different values oi Uf as a function of 
the c-c interaction Uc among conduction electrons. 

terms of the renormalization factors as 

A = -Zc + ^Zc^+AZcZfV^, (14) 

where Vg = V + g \s the effective hybridization defined 
before. It is seen from this formula that both of Zf 
and Zc have a tendency to reduce the gap A. For ex- 
ample, when Uc = 0, the gap is decreased by [//, lead- 
ing to the Kondo insulator with a small gap ~ AZfV"^. 
As mentioned above, however, the interplay of Uc and 
Uf increases the gap. We find that this enlargement 
is due to the c-f non-diagonal self-energy. Namely, in 
contrast to the simple PAM, the hybridization gap is af- 
fected not only hy Zf and Zc but also by the shift in 
the c-f hybridization, g, via the relation Vg = V -\- g vn 
the formula (14). Note that g is finite only in the case 
Uc and Uf ^ 0, and gives rise to the large gap 
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in Fig. 3. In this way, the Kondo insulator is adiabati- 

cally connected to the Mott insulator realized in the case 
Uc,Uf >> t, where the gap is mainly determined by [// 
or Uc- These results are consistent with those obtained 
by non-crossing approximation.^^) 

§4. Two-Orbital Periodic Anderson Model 

We now move to the two-orbital PAM to clarify the 
role of the orbital degrees of freedom. 

4..1 Correlations due to f-f interactions 

We first turn off the interactions for conduction elec- 
trons, Uc = U'^ = Jc = 0, and explore how the interac- 
tions for / electrons affect the formation of the Kondo 
insulator. Shown in Fig. 4 is the DOS for / electrons 



1 J 


. J 








uv=o.c 

I. 1 


...iliiJ 


ll J 


A 






UV=2.^ 

. XjX 111. 


. 1 i 1 




A 




1 ,11 


U'f=5.C 
111. 


,1 1 l.L 


, 1 1 liil 1 


111 


J Lil 11, 


UV=6.5 
111, 



-5 5 

CO 



Fig. 4. DOS of / electrons for different values of the inter-orbital 
/-/ Coulomb repulsion V'^ with a fixed Uf = 5.0. 

when the inter-orbital /-/ interaction U^- is varied with 
Uf = 5.0 being fixed. Increasing Up the upper and 
lower Hubbard- type bands, which are represented by 6- 
function like peaks, arc shifted to the higher energy re- 
gion with their weight being increased. On the other 
hand, qualitatively different behavior emerges for the gap 
formation in the lower energy region. Namely, the size 
of the gap is once enhanced with the increase ofU'j:, but 
beyond C/j ~ Uf it starts to decrease. This nonmono- 
tonic behavior is explicitly observed in Figs. 5 (a) and 
(b), where the renormalization factor Zf and the size of 
the gap A are plotted as a function of C/^. In the case 
Uf = 0, Zf and A decrease monotonically with the in- 
crease of U'j. The resulting insulating phase is regarded 
as a variant of the Kondo insulator, for which enhanced 
orbital (instead of spin) fluctuations reduce the gap. On 
the other hand, as already mentioned, a finite U / leads 
to the maximum structure in Zf and A around Uf ~ Uf. 
Therefore, low-energy electronic properties are quite sen- 
sitive to the balance of the Coulomb interactions U / and 
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Fig. 5. (a) R,cnormalization factor Zf, (b) spectral gap A, (c) 
and (d) local correlation functions for / electrons, as a function 
of the inter-orbital /-/ interaction Ui 



Up as pointed out in a related context of the Hubbard 
model. To confirm this, we also calculate the lo- 
cal correlation functions shown in Figs. 5 (c) and (d). 
The inter-orbital /-/ interaction Uj increases the double- 
occupation probability {n{^^^n^^^^) in the same orbital , 
while it decreases {nlg^nlf^} for two electrons occupying 
different orbitals. It should be noticed that these quanti- 
ties are altered dramatically around Uf ~ U'j.-, implying 
that orbital fluctuations are enhanced there. Therefore, 
we can say that the Kondo-insulating gap is enlarged due 
to orbital fluctuations among / electrons. 

We next consider the effects of the exchange coupling 
Jf to further clarify the role of local spin and orbital de- 
grees of freedom. In Figs. 6 (a) and (b), we show Zf 
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Fig. 6. (a) Renormalization factor Zf, (b) spectral gap A, and 
(c) local correlation functions as a function of the Hund coupling 
Jf. Other parameters are chosen as Uf = 10.0, i7c = = 0.0 
and Jc = 0.0. 

and A, which exhibit analogous J/ dependence. Since 
the Hund coupling restricts the available phase space of 
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/ electrons at each site, it has a tendency to suppress 
orbital fluctuations, thereby reducing both of Zf and A. 
Since the effective internal degrees of freedom depend on 
the sign of Jy, the above effects appear differently in two 
cases. When J j > 0, the Hund coupling favors a triplet 
state at each site, while the negative coupling J/(< 0) 
a singlet state. Therefore, for Jj < 0, orbital fluctua- 
tions are suppressed rather strongly, thus leading to a 
more prominent decrease in Zf and A. If U'j.- gets large, 
the system favors the local configuration of two electrons 
occupying the same orbital, so that the gap becomes in- 
sensitive to the Hund coupling in that parameter region, 
as seen in Figs. 6 (a) and (b). 

For reference, we show the local correlation functions 
(n{|n2|) and {n^T^n^^) in Fig. 6 (c), from which we can 
check to what extent the exchange coupling Jf affects 
spin and orbital fluctuations. It is seen that these quan- 
tities characterize the enhancement of triplet or singlet 
correlations depending on the sign of J/ , but only in the 
region Uj < Uf. In the case Uj > Uf, however, the dou- 
ble occupancy of either of two orbitals is favored, and 
thus {n{^^n(^^) gets larger than {n{.^n2-^) and {n{^n2^). 
This implies that the Hund coupling J/ is irrelevant in 
this case, consistent with the above results. 

4-2 Correlations due to c-c interactions 




Fig. 7. (a) Renormalization factor Zc, (b) eflfective hybridization 
Vg and (c) spectral gap A as a function of the intra-orbital c-c 
interaction Uc for different values of the inter-orbital /-/ inter- 
action U'jy. The intra-orbital /-/ Coulomb interaction is fixed as 
Uf = 5.0. 



We now study the influence coming from correlations 
among conduction electrons. In Figs. 7 and 8, several 
quantities discussed so far are plotted in the presence of 
the c-c Coulomb interaction Uc together with /-/interac- 
tions Uf and Uf. We notice that the effect of Uc shows 
up differently in the cases of Uf > Uf and Uf < Uf. 
For Uf > Uf, the gap A increases, although the renor- 
malization factor Zc decreases, reflecting electron corre- 
lations among conduction electrons (Figs. 7 (a) and (c)). 
This behavior is somewhat similar to the case discussed 
in the single-orbital model in Figs. 1 and 2. Namely, 
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Fig. 8. Local correlation functions as a function of intra-orbital 
c-c interaction Uc. The parameters are chosen as in 7 



the increase of the gap is caused by the increase of the 
hybridization Vg = V + g due to the c-f self-energy. 
This is indeed confirmed in Fig. 7 (b). In this case, 
the system changes continuously from the Kondo insu- 
lator to the Mott insulator This is typically seen in the 
local correlation functions shown in Fig. 8 (a), where 
both of the double-occupation probabilities 
and {n'^-^n'^^) arc suppressed with the increase of Uc. 

On the other hand, in the case of Uf < Uf, the gap 
A decreases with the increase of Uc, which is caused by 
the decrease of Vg. As shown in Figs. 8 (d) and (e), 
the double-occupation probability of the same orbital at 
Uc = is enhanced not only for / electrons {{n!^-^n^^)) 
but also for conduction electrons {{nj^-^n'^^)), even if 
C/^ = 0. This is a nontrivial effect due to the interplay of 
these interactions, and implies that the system atUc = 
is a variant of the Kondo insulator, for which the strong 
renormalization is caused by orbital fluctuations due to 
the inter-orbital interactions. The increase of Uc in turn 
suppresses such enhanced orbital fluctuations, as typi- 
cally seen in Figs. 8 (d) and (e), and at the same time 
gives rise to the decrease of Vg. Therefore, we may say 
that the resulting insulator is still in the Kondo insulat- 
ing phase with a small gap, which possesses both of the 
enhanced spin and orbital fluctuations. 

We finally observe what happens when we further 
introduce the inter-orbital c-c Coulomb interaction Uc 
among conduction electrons. We fix the values Uf = 5.0 
and Uc = 2.0, and compute the gap and the renormaliza- 
tion factor of conduction electrons as a function of U'. for 
various choices of Uf. Since we can follow the arguments 
outlined above, we only mention some interesting points 
briefly. As seen from Fig. 9, the gap is not so sensitive 
to the value of Uc for small Uf. A remarkable point is 
that when Uf is large (e.g. Uf = 6.0), all of Zc, Vg and 
A change dramatically as a function of U'c. Although the 
spectral gap A increases similarly to the case discussed in 
Fig. 7 for small Uf, we should notice that the nature of 
the insulating phase is different from each other. The dif- 
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Fig. 9. (a) Renormalization factor of conduction electrons Zc, (b) 
effective hybridization Vg , (c) spectral gap A as a function of the 
inter-orbital c-c Coulomb interaction V'^. We set Uf = 5.0 and 
Uc = 2.0 for simplicity. 
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Fig. 10. Local correlation functions for / and conduction elec- 
trons. We set Uf = 5.0 and Uc = 2.0. 



ference is distinguished in terms of the local correlation 
functions shown in Fig. 10. In particular, from panel (d) 
for Uj — 6.0, we can see that {n^-^n^^) and {n^-^n'^^) 

are enhanced, whereas (^{^nj^) and (TiiaTi2/3) ^^P' 
pressed. This implies that the insulating phase in that 
parameter region may be regarded as a variant of the 
Mott insulator with a large gap, for which orbital fluc- 
tuations arc enhanced while spin fluctuations are sup- 
pressed. This phase is qualitatively different from the 
one discussed in Fig. 7 for small U'j^, where spin (orbital) 
fluctuations are enhanced (suppressed) in the presence 

of C/e. 

Before closing this section, brief comments are in or- 
der for the impacts due to other interactions we have not 
addressed here. We have refrained from discussing the 
effects due to the exchange coupling among conduc- 
tion electrons, since they arc mostly the same as those 
of Jf. Also, we have not argued the influence due to the 
c-/ Coulomb repulsion in order to avoid the model to be 
too complicated. We have checked that this interaction 



has a tendency to enlarge the spectral gap. 

§5. Summary 

We have studied electron correlations for the orbitally 
degenerate periodic Anderson model by means of a lin- 
earized version of dynamical mean field theory. In par- 
ticular, we have focused on the role played by various 
intra- and inter-orbital interactions. 

By taking the single orbital model, we have first 
demonstrated how the c-c Coulomb repulsion naturally 
interpolates the Kondo insulator and the Mott insulator. 
In this case, the c-f part of the self-energy plays a ma- 
jor role rather than the renormalization effects due to 
and Zf. 

In the two orbital model, there are some remarkable ef- 
fects, which do not appear in the single-orbital case, due 
to the interplay of the intra- and inter-orbital interac- 
tions. One of the most interesting results is that orbital 
fluctuations are enhanced around Uf = Up which in 
turn suppress the renormalization effect. As a result, the 
Kondo insulating gap shows a maximum around there. 
This effect is somehow obscured by the exchange cou- 
pling Jf between / orbitals since it has a tendency to re- 
duce orbital fluctuations. Accordingly, the introduction 
of Jf again enhances the renormalization effect, making 
the Kondo gap smaller. 

Upon introducing the interactions for conduction elec- 
trons as well as / electrons in the degenerate model, there 
appear a variety of remarkable effects due to the interplay 
of these interactions. In fact, reflecting the subtle bal- 
ance of the interactions, the insulating phase with strong 
correlations has either the enlarged gap or the reduced 
gap, for which either of spin or orbital fluctuations (both 
of them in some cases) are enhanced. 

In this paper, we have restricted our analysis to the 
non-magnetic phase. We have indeed observed a num- 
ber of interesting properties even in such nonmagnetic 
insulating phase. It remains an interesting problem in 
the future study to clarify the competition between the 
nonmagnetic insulating state and magnetically ordered 
states, which may be particularly important in the sys- 
tem possessing large interactions. 

Acknowledgements 

The authors thank T. Mutou for valuable discussions. 

This work was partly supported by a Grant-in- Aid from 
the Ministry of Education, Science, Sports and Culture 
of Japan. A part of computations was done at the Super- 
computer Center at the Institute for Solid State Physics. 
University of Tokyo and Yukawa Institute Computer Fa- 
cility. 



1) T. Susaki et at Phys. Rev. Lett. 77 (1996) 4269. 

2) K. Sugiyama, F Iga, and M. Kasaya: J. Phys. Soc. Jpn. 57 
(1988) 3946. 

3) T. Saso and H. Harima: J. Phys. Soc. Jpn. 68 (1999) 2491. 

4) J. W. Allen, B. Batlogg, and P. Wachter: Phys. Rev. B. 20 

(1979) 4807. 

5) M. F. Hundley, P. C. Canficld, and Z. Fisk: Phys. Rev. B. 42 
(1990) 6842. 

6) A. Severing, J. D. Thompson, P. C. Canfield, Z. Fisk, and P. 



Periodic Anderson model with degenerate orbitals 



7 



S. Riscboiough: Phys. Rev. B 44 (1991) 6832. 

T. E. Mason et al: Phys. Rev. Lett. 69 (1992) 490. 

P. S. Riseborough: Phys. Rev. B 45 (1992) 13984. 

C. Sanchez-Castro, K. S. Bedell, and B. R. Cooper: Phys. 
Rev. B 47 (1993) 6879. 

J. C. Coolcy, M. C. Aronson, Z. Fisk, and P. C. Canfield: 
Phys. Rev. Lett. 74 (1995) 1629. 

T. Brugger, T. Schrcincr, G. Roth, P. Adelmann, and G. 

Czjzck: Phys. Rev. Lett. 71 (1993) 2481. 

Y. M. Li: Phys. Rev. B. 52 (1995) 6979. 

J. Igarashi, K. Murayama, and P. Pulde: Phys. Rev. B. 52 

(1995) 15966. 

K. Itai and P. Fazekas: Phys. Rev. B. 54 (1996) 752. 

T. Schork and S. Blawid: Phys. Rev. B. 56 (1997) 6559. 

Y. One, T. Matsuura, and Y. Kuroda:J. Phys. Soc. Jpn. 63 

(1993) 1406. 

M. Jarrell, H. Pang, D. L. Cox, F. Anders, and A. Chattopad- 
hyay: Physica B 230-232 (1997) 557. 

A. Tsuruta, A. Kobayashi, and Y. Ono: J. Phys. Soc. Jpn. 
68 (1999) 2491. 

W. Metzner and D. VoUhardt: Phys. Rev. Lett. 69 (1989) 
324. 

E. Miiller-Hartmann, Z. Phys. B: Condens. Matter 74 (1989) 
507. 

T. Pruschke, M. Jarrell, and J. K. Freericks: Adv. Phy. 42 
(1995) 187. 

A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg: 
Rev. Mod. Phys 68 (1996) 13. 

A. Georges and G. Kothar: Phys. Rev. B. 45 (1992) 6479. 
M. Jarrell: Phys. Rev. Lett. 69 (1992) 168. 
X. Y. Zhang, M. J. Rozenberg, and G. Kotliar: Phys. Rev. 
Lett. 70 (1993) 1666. 

M. Caffarel and W. Krauth: Phys. Rev. Lett. 72 (1994) 1545. 

D. S. Fisher, G. Kothar, and G. Moeller: Phys. Rev. B. 52 

(1995) 17112. 

H. Kajueter and G. Kotliar: Phys. Rev. Lett. 77 (1996) 131. 

R. Bulla: Phys. Rev. Lett. 83 (1999) 136. 

M. Potthoff: Phys. Rev. B. 64 (2000) 165114. 

R. Bulla, and M. PotthofT: Eur. Phys. J. B. 13 (2000) 257. 

Y. Ono, R. Bulla, A. C. Hewson, and M. PotthofT: Eur. Phys. 

J. B. 22 (2001) 283. 

M. J. Rozenberg: Phys. Rev. B. 55 (1997) R4855. 

K. Held and D. VoUhardt: Euro. Phys. J. B 5 (1998) 473. 

J. E. Han, M. Jarrell, and D. L. Cox: Phys. Rev. B. 58 (1998) 

R4199. 

Th. Maier, M. B. Zolfl, Th. Pruschke, and J. Keller: Eur. 
Phys. J. B 19 (1999) 377. 

T. Momoi and K. Kubo: Phys. Rev. B. 58 (2000) R567. 
Y. Imai and N. Kawakami: J. Phys. Soc. Jpn. 70 (2001) 2365. 
A. Koga, Y. Imai, and N. Kawakami: Phys. Rev. B. 66 (2002) 
165107. 

V. S. Oudovenko and G. Kotliar: Phys. Rev. B. 65 (2002) 
0750102. 

S. Florens, A. Georges, G. Kotliar, and O. ParcoUet: Phys. 
Rev. B. 66 (2002) 205102. 

Y. Ono, M. Potthoff, and R. Bulla: Phys. Rev. B. 67 (2003) 
035119. 

A. Koga, T. Ohashi, Y. Imai, S. Suga and N. Kawakami: J. 

Phys. Soc. Jpn. 72 (2003) 1306. 

A. Liebsch: Phys. Rev. Lett. 91 (2003) 226401. 

A. Koga, N. Kawakami, T.M. Rice, and M. Sigrist: cond- 

mat/0401223. 

T. Mutou and D. Hirashima: J. Phys. Soc. Jpn. 63 (1994) 

4475; T. Mutou: Phys. Rev. B. 64 (2001) 165103. 

M. Jarrell: Phys. Rev. B. 51 (1995) 7429. 

T. Saso and M. Itoh: Phys. Rev. B. 53 (1996) 6877. 

Th. Pruschke, R. Bulla, M. Jarrell: Phys. Rev. B. 61 (2000) 

12799. 

D. Meyer and W. Nolting: Phys. Rev. B. 61 (2000) 13465. 
Chang-Il Kim, Y. Kuramoto, and T. Kasuya: J. Phys. Soc. 
Jpn. 59 (1990) 2414. 



